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(54) Lautheitsgesteuerte Verarbeitung akustischer Signale 



(57) Mit dem Verfahren werden akustische Signale 
bspw. in Horgeraten auf eine solche Weise lautheitsge- 
steuert verarbeitet, dass die vom Horbehinderten sub- 
jektiv empfundene Lautheit wieder stets der von Nor- 
malhorenden empfundenen Lautheit entspricht. Die Si- 
gnatverarbeitung erfolgt ohne Fouriertransformation 
und ohne Aufteilung des Signals in Teilbandsignale auf 
iterative Weise vollstandig im Zeitbereich. Damit wird 
der Nachteil inakzeptabel langer Signalverzogerungs- 
zeiten bekannter Verfahren eliminiert und eine prakti- 



sche Nutzung Oberhaupt erst ermoglicht. Die Vorrich- 
tung zur DurchfOhrung des Verfahrens enthalt eine Ver- 
arbeitungsstufe (4) zur iterativen Berechnung einer fur 
die Lautheit charakteristischen Steuergrdsse (y) und ei- 
ne damit zeitabhangig gesteuerte Korrekturfilterstufe 
(7). Das erfindungsgemasse Verfahren erfordert gegen- 
Ober bekannten Verfahren nurdrastisch reduzierte Ver- 
arbeitungsressourcen, was hauptsachlich auf die be- 
sonders effiziente und unkonventionelle Ausfuhrung der 
Verarbeitungsstufen zuruckzufuhren ist. 
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Beschreibung 

Die Erfindung betrifft ein Verfahren zur lautheitsgesteuerten Verarbeitung akustischer Signale in Schallverarbei- 
tungsgeraten sowie eine Vorrichtung zur Durchfuhrung des Verfahrens gemass den Oberbegriffen der unabhangigen 

5 Patentanspruche. Die Erfindung eignet sich besonders fur den Einsatz in Horgeraten fur Hdrbehinderte; eintreffende 
akustische Signale werden auf eine solche Weise verarbeitet, dass die vom Horbehinderten subjekliv empfundene 
Lautheit wieder stets der von Normalhorenden empfundenen Lautheit entspricht. 

Die Idee der lautheitsgesteuerten Verarbeitung akustischer Signale ist sett langerem bekannt und von verschie- 
denen Autoren in der Fachliteratur beschrieben worden, so z.B. auch von N. Dillier et al. im "Journal of Rehabilitation 

10 Research and Development", Vol. 30, No. 1, 1993, S. 100-103. Das Verfahren beruht darauf, dass Normalhorenden 
und Horbehinderten bekannte Testsignale zur Beurteilung der subjektiv empfundenen Lautheit vorgefuhrt werden. Als 
Testsignale werden harmonische Sinussignale Oder Schmalbandrauschen verwendet. Die subjektiv empfundene Laut- 
heit hangt von der Signalleistung und von der Frequenz eines Sinussignals, bzw. von der Frequenz der dominanten 
Signalkomponenten eines komplexen Signals, ab. Die subjektiven Lautheitsangaben werden auf einer normierten 

'5 Skata mit Wertebereich [0, 1] erfasst. Durch Vergleich der Angaben eines Horbehinderten mit jenen einer Referenz- 
gruppe Normalhorender lassen sich horbehindertenspezifische, lautheitsabhangige Korrekturdaten ermitteln. In einem 
dazu passenden Signalverarbeitungsverfahren werden diese Korrekturdaten dann eingesetzt, um dem Horbehinderten 
die akustischen Signale seiner Umwelt in der zielsetzungsgemassen Art aufzubereiten. Im erwahnten Artikel wurden 
in Verstandlichkeitstests mit einer Gruppe von 13 Horbehinderten bemerkenswerte Verbesserungen der Verstandlich- 

20 keit nachgewiesen. 

Trotz der audiologischen Wirkung kann die lautheitsgesteuerte Verarbeitung in der bisher bekannten Form in der 
Praxis nicht eingesetzt werden. Wie im erwahnten Artikel beschrieben, erfolgt namlich die Verarbeitung durch Fourier- 
transformation kurzer Signalsegmente, Modifikation der Kurzzeitspektren und Rucktransformation der modifizierten 
Kurzzeitspektren in den Zeitbereich. Infolge der segmentweisen Verarbeitung ergibt sich fur das verarbeitete Signal 

25 eine Verzogerung von nahezu 20 ms. Bei Verstandlichkeitstests spielt diese Verzogerung keine Rolle. In der Praxis, 
wenn der Horbehinderte selbst auch spricht und seine eigene Stimme dermassen verzogert wahrnimmt, ist sie jedoch 
vollig inakzeptabel. Bei dem im erwahnten Artikel beschriebenen Verfahren betragt die Dauer der einzelnen Segmente 
12.8 ms, und dieser Wert kann auch nicht wesentlich unterschritten werden, weil zur Gewinnung eines brauchbaren 
Kurzzeitspektrums eine minimale Segmentdauer in dieser Grossenordnung unerlasslich ist. 

30 Als Alternative zur segmentweisen Verarbeitung wurde gelegentlich der Ansatz versucht, das akustische Signal 

in Teilbandsignale aufzuteilen und die einzelnen Teilbandsignale mit separaten Verstarkungswerten zu verarbeiten. 
Aus praktischen Erprobungen ist bekannt, dass sich beim Autteilen in bis zu drei Teilbandsignale Verbesserungen 
erzielen lassen. Eine Aufteilung in mehr Teilbandsignale fOhrt jedoch wieder zu schlechteren Ergebnissen. Ein Grund 
dafur sind moglicherweise die an den Teilbandgrenzen entstehenden Diskontinuitaten der Ubertragungsfunktion. Stellt 

35 man der Aufteilung des Signals in drei Teilbandsignale die Frequenzauflosung der Kurzzeitspektren der segmentweisen 
Verarbeitung gegenuber, so ist auf jeden Fall klar, dass das Potential der segmentweisen Verarbeitung mit dem alter- 
nation Ansatz nicht ausgeschopft werden kann. Und selbst wenn sich mit dem Aufteilen in mehr Teilbandsignale Wege 
zu verbesserten Ergebnissen finden liessen, entstOnde dabei wiederum das Problem deutlich zunehmender Signal- 
verzogerung. 

40 Ein weiterer Aspekt fur das gute Gelingen der lautheitsgesteuerten Signalverarbeitung hat mit dem in der Verar- 

beitung verwendeten Lautheitsmodell zu tun. Im Gegensatz zu den einfachen Testsignalen verteilt sich die Signallei- 
stung von Sprache, Musik und Gerauschen zeitabhangig und in komplexer Art uber ein weites Frequenzintervall. Mit 
einem Lautheitsmodell wird diesen komplexen Signalen zeitabhangig ein Lautheitswert zugeordnet, der im Ideatfall 
exakt mit der von Normalhorenden empfundenen Lautheit zusammenfallt. Der mit dem Lautheitsmodell ermittelte Wert 

45 wird fur die zeitabhangige Steuerung der Signalverarbeitung eingesetzt. Das im erwahnten Artikel beschriebene Laut- 
heitsmodell berucksichtigt nebst der Gesamtenergie eines Signalsegments auch die Schwerpunktsfrequenz seines 
Kurzzeitspektrums. Zur Berechnung der Schwerpunktsfrequenz wird auf Grundlagen von E. Zwicker zuruckgegrrffen, 
die in dessen Lehrbuch "Psychoakustik", Springer Verlag, Berlin, 1982, auf den Seiten 51 bis 53 zusammengefasst 
sind. Aus den Spektrallinien des Kurzzeitspektrums werden in einem ersten Schritt die Energien E(z) der einzelnen 

50 Frequenzgruppen gebildet und anschliessend in Analogic zur Berechnung des Schwerpunkts in der Mechanik zu einer 
Schwerpunktsfrequenz 

c = Zz-E(z)/2:E(z) (1) 

55 

auf der mit z bezeichneten Barkskala verrechnet. Wollte man dieses Lautheitsmodell mrttels Aufteilung des Signals in 
Teilbandsignale realisieren, so waren zur Verarbeitung einer Bandbreite von 7700 Hz insgesamt 21 Teilbandsignale 
unterschiedlicher Bandbreite entsprechend den bekannten Frequenzgruppenbreiten zu bilden. Nebst der bereits er- 
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wahnten, stark ansteigenden Signalverzogerung wurde dieses Vorgehen auch ausserordentlich grosse rechnerische 
Ressourcen erfordern. Mit den derzeit verfugbaren Technologien fur integrierte Schaltungen ist damit aber, wie auch 
fur den Ansatz mit segmentweiser Verarbeitung, die Umsetzung in ein Horgerat mit heute ublichen geometrischen 
Abmessungen und Stromverbrauch ausgeschlossen. 

5 Es ist Aufgabe der vorliegenden Erfindung, ein Verfahren zur lautheitsgesteuerten Verarbeitung akustischer Si- 

gnale in Schallverarbeitungsgeraten anzugeben, welches insbesondere in Horgeraten einsetzbar ist. Die vom Horge- 
ratbenutzer subjektiv empfundene Lautheit soli stets der von Normalhorenden empfundenen Lautheit entsprechen. 
Insbesondere soil die Signalverzogerung so klein sein, dass ein Horgeratbenutzer durch die verzogerte v\fahmehmung 
seiner eigenen Stimme beim Sprechen nicht irritiert wird. Es sollen auch die rechnerischen Ressourcen gegenuber 

10 bekannten Verfahren fur die lautheitsgesteuerte Verarbeitung akustischer Signale reduziert werden. Weiter soil eine 
Vorrichtung zur Durchfuhrung des erfindungsgemassen Verfahrens geschaffen werden. 

Die Aufgabe wird gelost durch das Verfahren und die Vorrichtung gemass den unabhangigen Patentanspruchen. 
Im erfindungsgemassen Verfahren erfolgt die Verarbeitung des akustischen Signals ohne Fouriertransformation, 
also vollstandig im Zeitbereich, und auch ohne Aufteilung in Teilbandsignale. Die Besonderheit des erfindungsgemas- 

15 sen Verfahrens liegt darin, dass eine fur die Lautheit charakteristische Steuergrosse \\t auf iterative Weise berechnet 
und zur Steuerung eines zeitabhangigen Korrekturfi Iters eingesetzt wird. Mit dem Ausdruck "iterative Berechnungs- 
weise" ist gemeint, dass fur die Steuergrosse \j/ zu jedem Abtastzeitpunkt ein neuer Wert berechnet wird, und zwar 
unter Verwendung von Werten. welche die zu ihrer Berechnung notwendigen Grossen im jeweils vorhergehenden 
Abtastzeitpunkt aufwiesen. Im Gegensatz zum bekannten segmentweisen Verfahren wird die lautheitsspezifische 

20 Steuergrosse also nicht nur als Mittetwert aufeinanderfolgender Signalsegmente, sondern als kontinuierliche Zeitfunk- 
tion ermittelt. Die kurze, typischerweise auf 2 ms bemessene Signalverzogerung stellt die fur eine zuverlassige Schatz- 
wertbildung uber den jeweiligen Gultigkeitszeitpunkt hinaus erforderliche Beobachtungsdauer dar und ist somit, im 
Gegensatz zum segmentweisen Verfahren, nicht bloss die Folge einer nachteiligen Eigenschaft der gewahlten Imple- 
mentation. Die iterative Berechnungsweise erfolgt im erfindungsgemassen Verfahren mittels besonders effizienter und 

25 zugleich origineller Verfahrensschritte. 

Das zeitabhangige Korrekturfi Iter wird dadurch gesteuert, dass Parametern des Korrekturfi Iters zu jedem Abtast- 
zeitpunkt durch Interpolation mit Hilfe der Steuergrosse y neue Werte zugewiesen werden. Im Gegensatz zum seg- 
mentweisen Verfahren, wo die horbehindertenspezifischen Korrekturdaten als Verstarkungswerte fur die einzelnen 
Spektrallinien eines Kurzzeitspektrums gespeichert sind, werden beim erfindungsgemassen Verfahren fur wohldefi- 

30 nierte Werte der Steuergrosse \j/ Koeffidentensatze fur Prototypenfilter im voraus bestimmt und gespeichert. Dabei 
verlaufen die Ubertragungsfunktionen dieser Prototypenfilter entlang den entsprechenden Verstarkungswerten, die 
beim segmentweisen Verfahren fur die einzelnen Spektrallinien eines Kurzzeitspektrums bestimmt sind. Zur Charak- 
terisierung der Prototypenfilter werden im erfindungsgemassen Verfahren Koeffizientensatze verwendet, von denen 
bekannt ist, dass sie sich fur eine Interpolation eignen, d.h. dass die durch interpolierte Koeffizienten bestimmte Uber- 

35 tragungsf unktion erwartungsgemass zwischen den Ubertragungsfunktionen verlauft, welche durch die der Interpolation 
zugrunde liegenden Koeffizientensatze bestimmt sind. 

Mit dem erfindungsgemassen Verfahren werden also vollig neue Wege beschritten. Es werden die im erwahnten 
Artikel von N. Dillier et al. beschriebenen guten Verstandlichkeitsergebnisse erzielt. Daruber hinaus reduziert aber das 
erfindungsgemasse Verfahren die Signalverzogerung auf ca. 2 ms und erreicht zugleich eine drastische Reduktion 

40 der rechnerischen Ressourcen. Daher ist es moglich, das erfndungsgemasse Verfahren in einem Horgerat heute ub- 
licher Bauform umzusetzen. 

Die Erfindung betrifft weiter eine Vorrichtung zur Durchfuhrung des erfindungsgemassen Verfahrens. Diese Vor- 
richtung enthalt eine Stufe zur iterativen Berechnung der fur die Lautheit charakteristischen Steuergrosse y und eine 
damit zeitabhangig gesteuerte Korrekturfilterstufe, welche eintreffende akustische Signale zielsetzungsgemass auf- 
45 bereitet. Die bereits erwahnte drastische Reduktion der benotigten Verarbeitungsressourcen hat verschiedene Ursa- 
chen. Zum einen entfallt bei der iterativen Berechnungsweise die segmentweise Pufferung des Eingangs- und Aus- 
gangssignals. Dann ergibt sich beim Speichem der Koeffizientensatze fur die Prototypenfilter auch eine wesentliche 
Einsparung gegenuber dem Speichern der Verstarkungswerte fur die einzelnen Spektrallinien der Kurzzeitspektren. 
Im folgenden wird die Erfindung anhand der Zeichnungen und eines detaillierten Ausfuhrungsbeispiels naher er- 
50 lautert. Es zeigen: 



Fig. 


1 


ein Blockdiagramm der lautheitsgesteuerten Verarbeitung im Uberblick, 


Fig- 


2 


ein Blockdiagramm zur Ermittlung der fur die Lautheit charakteristischen Steuergrosse, 


Fig. 


3 


ein Signalflussdiagramm eines rekursiven Digitalfilters, 


Fig. 


4 


ein Signalflussdiagramm einer einfachen Schatzwertberechnungseinheit, 
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20 



Fig. 5 ein Signalflussdiagramm einer Schatzwertberechnungseinheit fur die Signalleistung, 

Fig. 6 und 7 Schemas zur Gewinnung von Tabellenadressen, 

Fig. 8 ein Signalflussdiagramm einer Schatzwertberechnungseinheit fur die Schwerpunktsfrequenz, 

Fig. 9 ein Signalflussdiagramm eines nichtlinearen Glattungsfiiters, 

Fig. 10 ein Diagramm fur den Zusammenhang interner Grossen des nichtlinearen Glattungsfiiters, 

Fig. 11 ein Schema fur eine zweidimensionale Interpolation, 

Fig. 12 ein Blockdiagramm der Interpolation von Parametern des Korrekturf liters, 

15 Fig. 13 ein Schema zur Gewinnung von Tabellenadressen und Proportionalgrossen fur Interpolationen, 

Fig. 14 ein Blockdiagramm des zeitabhangigen Korrekturfilters, 

Fig. 15 ein Signalflussdiagramm eines Kreuzgliedfilters zur Realisierung von Nullstellen, 

Fig. 16 ein Signalflussdiagramm eines Kreuzgliedfilters zur Realisierung von Polstellen, 

Fig. 17 und 18 Schemas fur zweistufige lineare Interpolationen und 

25 Fig. 19 und 20 Schemas zur Gewinnung von Tabellenadressen und Proportionalgrossen fur Interpolationen. 

Figur 1 zeigt den Einsatz des erfindungsgemassen Verfahrens und das Verfahren selbst in einem schematischen 
Uberblick. Ein akustisches Signal wird von einem Mikrofon 1 in ein elektrisches Signal gewandelt, das von einem 
Signalumsetzer 2 digitalisiert und anschliessend in einem Hochpassfilter 3 von allfalligem Offset und tiefstfrequenten 
30 Storsignaikomponenten befreit wird. 

Die wesentlichen Schritte des erfindungsgemassen Verfahrens bestehen in der Verarbeitung eines Ausgangssi- 
gnals x des Hochpassfilters 3. In einer Verarbeitungsstufe 4 erfolgt die iterative Berechnung der Steuergrosse vy. in 
einer anschliessenden Interpolationsstufe 5 werden damit die Parameter eines zeitabhangigen Korrekturfilters 7 er- 
mittelt und an dieses ubergeben. Eine Verzogerungsstufe 6 sorgt hinsichtlich der Filterung mit dem Korrekturfilter 7 
35 fur die Synchronisation des Signals x mit den aus ihm abgeleiteten Filterparameterwerten, indem sie eine entspre- 
chende Signalverzogerung, beispelsweise urn 2 ms, bewirkt. Bei einer Abtastrate von 16 kHz wird die Verzogerungs- 
stufe 6 vorteilhafterweise als zyklischer Puffer mit 32 Speicherplatzen ausgelegt. 

Das mit dem Korrekturfilter 7 gefilterte Signal y gelangt zu einem Signalumsetzer 8 und wird dort in ein analoges 
elektrisches Signal gewandelt. In einer analogen Verstarkerstufe 9 wird es noch mit einem horbehindertenspezifischen, 
40 aber zeitlich konstanten Verstarkungswert g e verstarkt und anschliessend einem elektro-akustischen Signalwandler 
10 zugefuhrt. Der Wert von g e wird beim Aufbereiten der Koeffizientensatze fur die Prototypenfilter bestimmt, und zwar 
so, dass das in der Vorrichtung zur Durchfuhrung des Verfahrens verwendete 16 Bit breite Zahlenformat moglichst 
optimal genutzt wird, wobei eine Begrenzung der verarbeiteten Signale infolge der in der Vorrichtung vorausgesetzten 
Sattigungsarithmetik jedoch nur ausnahmsweise wirksam werden soil. 
45 wie bereits erwahnt wurde, kann die Lautheit komplexer Signale aufgrund der Gesamtenergie kurzer Signalseg- 

mente und der Schwerpunktsfrequenz ihrer Kurzzeitspektren ermittelt werden. Dabei hangt die Lautheit naherungs- 
weise quadratisch von der auf einer logarithm ischen Skala ausgedruckten Signalenergie ab. Wie noch gezeigt wird, 
lasst sich das Lautheitsmodell im erfindungsgemassen Verfahren mit einer zweidimensionalen linearen Interpolation 
implementieren. Diese Interpolation liefert genauere Ergebnisse, wenn die Steuergrosse 



50 



V = NL'-VL min y(VL niax .VL m|n ) (2) 



eingefuhrt wird, die folglich naherungsweise linear von der logarithmischen Signalenergie abhangt. Dabei stellt U die 
55 auf den Wertebereich [L mln , L^] begrenzte Lautheit dar, und L mjn und L max sind sinnvoll gewahlte Minimal- und Ma- 
ximalwerte der Lautheit, die somit den Arbeitsbereich des Verfahrens definieren, innerhalb dessen das Korrekturfilter 
aufgrund kleinster Veranderungen der Lautheit stetig nachgefuhrt wird. Aufgrund der Formel (2) ist y eine auf den 
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Wertebereich [0, 1] normierte Steuergrdsse, und fur Lautheitswerte ausserhalb des Wertebereichs [L mjn , L ma J wird 
das Korrekturfilter fur y = 0 bzw. jenes fur y = 1 eingeselzt. 

Das Blockdiagramm in Fig. 2 zeigt etwas detaillierter, wie die Steuergrossey aus dem Eingangssignal x gewonnen 
wird. Verglichen mit dem bekannten segmentweisen Verfahren tritt beim erfindungsgemassen iterativen Signalverar- 
beitungsverfahren an die Stelle der Signalenergie eines kurzen Signalsegments eine momentane Signalleistung q und 
an die Stelle der Schwerpunktsfrequenz seines Kurzzeitspektrums eine momentane Schwerpunktsfrequenz c. Diese 
Grossen werden in Verarbeitungsstufen 11-15 ermittelt. Nach einer Verarbeitungsstufe 13 weisen entsprechende Aus- 
gangssignalwerte c r und q r aufgrund der iterativen Berechnungsart noch eine unerwunschte Streuung auf, die in nach- 
folgenden Glattungsfiftern 1 4 und 1 5 eliminiert wird. Die geglatteten Signale c und q werden in einer Verarbeitungsstufe 
16 der bereits erwahnten zweidimensionalen Interpolation zugefuhrt, wobei die aufeinanderfolgenden Ausgangssi- 
gnalwerte yr ebenfalls eine noch unerwunschte Streuung aufweisen, die mit einem nachfolgenden Glattungsfilter 17 
eliminiert wird. 

Ein wesentlicher Aspekt des erfindungsgemassen Verfahrens liegt in der iterativen Berechnungsart der logarith- 
mischen Signalleistung q sowie der auf einer Barkskala ausgedrOckten Schwerpunktsfrequenz c, also der Umsetzung 
der Formel (1) in ein iteratives Berechnungsschema. Anstelle der Bildung der f requenzgruppenspezifischen Energien 
E(z) wird im erfindungsgemassen Verfahren eine frequenzselektive Gewichtung des Eingangssignals x mit einem Filter 
vorgenommen, das im weiteren als Frequenzgruppenfilter bezeichnet wird. Das Frequenzgruppenfilter ist in Fig. 2 als 
Verarbeitungsstufe 1 1 dargestellt, und sein Ausgangssignal ist mit cp bezeichnet. Seine von der Frequenz f abhangige 
Ubertragungsfunktion 

H FG (f) = VAf Q (f)/VAf G (f N ) (3) 

ergibt sich aus der Frequenzgruppenbreitefunktion Af G (f). Der Nenner in Formel (3) bewirkt eine Normierung, wobei 
f N die Nyquistfrequenz, im Ausfuhrungsbeispiel also 8 kHz, bezeichnet. Die Normierung dient der optimalen Nutzung 
des im Ausfuhrungsbeispiel vorgegebenen, 16 Bit breiten Festkommazahlenformats. Im Ausfuhrungsbeispiel wird die 
Ubertragungsfunktion H FG (f) durch ein rekursives Filter 11 zweiter Ordnung approximiert. Die Struktur des Frequenz- 
gruppenfilters 11 ist der Vollstandigkeit halber in Fig. 3 illustriert 

Anstelle der Gewichtung der Frequenzgruppenenergien E(z) mit den Frequenzgruppenindizes z im Zahler von 
Formel (1 ) wird im erfindungsgemassen Verfahren eine frequenzselektive Gewichtung des Signals cp mit einem Filter 
vorgenommen, das im wertern als Barkfilter bezeichnet wird. Das Barkfilter ist in Fig. 2 als Verarbeitungsstufe 12 
dargestellt, und sein Ausgangssignal ist mit bezeichnet Seine Obertragungsfunktion 

H B (f) = Vz(f)/Vz(f N ) (4) 

ergibt sich aus der Tonheitsfunktion z(f). Der Nenner in Formel (4) bewirkt wiederum eine Normierung zum Zwecke 
einer optimalen Nutzung des vorgegebenen Zahlenformats. Im AusfOhrungsbeispiel ist die Ubertragungsfunktion H B 
(f) ebenfalls durch ein rekursives Dig'rtalfilter 12 zweiter Ordnung approximiert, das wiederum die in Fig. 3 abgebildete 
Struktur aufweist. 

Mit den Signalen d und cp ist es im erfindungsgemassen Verfahren mdglich, die momentane Schwerpunktf requenz 
gemass Formel (1) auf iterative Weise zu berechnen. Dazu wird in einer Verarbeitungsstufe 13 der Quotient ihrer 
Signalleistungen berechnet. 

Zur iterativen Berechnung von Signalleistungen wird im erfindungsgemassen Verfahren auf eine einfache Schatz- 
wertberechnungseinheit erster Ordnung fur den zeitlich exponentiell gewichteten Erwartungswert des quadrierten Ein- 
gangssignals zuruckgegriffen. Eine solche Schatzwertberechnungseinheit ist fur den allgemeinen Fall, mit Eingangs- 
signal u und Ausgangssignal v, in Fig. 4 dargestellt. In diesem Signalflussdiagramm ergibt sich ein neuer Ausgangs- 
signalwert v dadurch, dass der Ausgangssignalwert des vorhergehenden Abtastzeitpunkts mit der Konstanten (1 - e) 
multipliziert und zu diesem Produkt das mit dem konstanten Faktore multiplizierte Quadrat des neuen Eingangssignal- 
wertes u hinzuaddiert wird. Mit der Adaptionskonstanten e, fur welche 



0 < e < 1 

gilt, lasst sich die Geschwindigkeit steuern, mit der das Ausgangssignal v der sich verandernden Eingangssignallei- 
stung folgt. 
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Die einfache Schatzwertberechnungseinheit von Fig. 4 weist die Nachteile auf, dass zur Verarbeitung des qua- 
drierten Eingangssignals ein doppelt breites Zahlenformat erforderlich ist, und dass fur die nachfolgenden Berechnun- 
gen zusatzlich noch der Logarithmus des Ausgangssignals v benotigt wird. Beide Aspekte werden im erfindungsge- 
massen Verfahren auf einfache Art, wie in Fig. 5 gezeigt, durch Einbettung der einfachen Schatzwertberechnungsein- 

5 heit von Fig. 4 in eine digrtale Regelschleife gelost. 

Die Funktionsweise des Signalflussdiagramms in Fig. 5 beruht darauf, dass die Grosse v auf einen fest vorgege- 
benen Sollwert geregelt wird. Zu diesem Zweck wird fur jeden neu berechneten Signalwert v die inkrementelle log- 
arithmische Zuwachs- bzw. Abnahmegrosse der Signalleistung bestimmt, welche der Abweichung des Wertes v vom 
vorgegebenen Sollwert entspricht. Die gesuchte logarithmische Signalleistung p ergibt sich in der Folge durch blosses 

10 Akkumulieren der aufeinanderfolgenden inkrementellen Veranderungswerte. Fur das richtige Funktionieren der Re- 
gelschleife ist es notwendig, dass jeder Eingangssignalwert x mit einem jeweils zum Schatzwert p passenden Skalie- 
rungsfaktor skaliert wird und dass auch die Grosse v selbst vor dem erneuten Aufdatieren mit einem der Leistungs- 
anderung entsprechenden Justierwert multipiikativ nachgefuhrt wird. Die Ermittlung sowohl der inkrementellen Ande- 
rung wie auch der Skalierungs- und Justierwerte erfolgt im erfindungsgemassen Verfahren zu jedem Abtastzeitpunkt 

75 fQr Werte der Grossen v und p, deren Genauigkeit durch Abschneiden auf 6 bzw. 7 Nachkommastellen beschrankt ist. 
Damit ist der effiziente Einsatz von Tabellen moglich, in denen die 64 bzw. 128 im voraus berechneten, passenden 
Werte abgespeichert sind. Schliesslich brauchen zur Adressierung der Tabellen dann lediglich, wie in den Fig. 6 und 
7 gezeigt, die relevanten Bitfelder aus den Grossen v und p extrahiert zu werden. In Fig. 5 ist die Tabelle mit den 
inkrementellen logarithmischen Leistungsanderungen mit Ap bezeichnet Zum Einsparen ansonsten separat auszufuh- 

20 render Multiplikationen enthalt die Tabelle S in Fig. 5 ferner noch modifizierte Skalierungswerte, die aus den ursprung- 
lichen Skalierungswerten durch Multiplikation mit der Wurzel aus der Konstanten e gewonnen wurden. Zum gleichen 
Zwecke sind die Justierwerte in der mit A bezeichneten Tabelle bereits mit der Konstanten (1 - e) multipliziert. Ferner 
reicht fur die Speicherung der Grossen v und p sowie aller Tabellenwerte in Fig. 5 das ubliche, 16 Bit breite Festkom- 
mazahlenformat. 

25 wie bereits erwahnt wurde, beruht im erfindungsgemassen Verfahren die iterative Berechnung der Schwerpunkts- 

frequenz auf der Berechnung des Quotienten der Signalleistungen der Signale O und cp, beispielsweise in der Verar- 
beitungsstufe 13. Die Berechnung der Signalleistungen wird auf das in Fig. 5 abgebildete Signatflussdiagramm zu- 
ruckgef uhrt. Damit ergibt sich fur die Berechnung der Schwerpunktsf requenz das in Fig. 8 dargestellte Signalflussdia- 
gramm. Der untere Teil des Diagramms ist identisch mit der Fig. 5. Er dient der Berechnung der Leistung des Signals 

30 cp. Der obere Teil dient der Berechnung der Leistung des Signals Bei dieser Berechnung werden die Skalierungs- 
und Justierwerte aus dem unteren Schaltungsteil ubernommen, wodurch sich das Signatflussdiagramm im oberen Teil 
gegenuber Fig. 5 vereinfacht. Mit dieser Anordnung ist die optimale Nutzung des Zahlenformats auch fOr die Berech- 
nung der Leistung des Signals 0 gewahrleistet, und die gesuchte Schwerpunktsfrequenz ergibt sich, wie erwahnt, 
durch Quotientenbildung der beiden Signalleistungen. 

35 Wie in Fig. 8 gezeigt, erfolgt die Berechnung eines aus einem Zahler Z und einem Nenner N gebildeteten Quoti- 

enten Q = Z/N anhand der bereits mit einem Justierwert aus der Tabelle A nachgefuhrten Signalleistungswerten. Das 
hat den Vorteil, dass die ansonsten benotigte, ungunstige Division signifikant vereinfacht werden kann. In einem auf 
den vorgegebenen Sollwert normierten Zahlenformat nimmt der Nenner 

40 

N = 1 + 6 mit |6| <* I (6) 

nur unwesentlich von 1 verschiedene Werte an, und anstelle der Division durch (1 + 8) kann der Quotient 

45 

Q - Z.<1 - 8) (7) 

durch Multiplikation des Zahlers Z mit (1 - 8) angenahert werden. 

50 Wie schon zuvor ausgefuhrt wurde, kann die Lautheit aus der Signalleistung p und der Schwerpunktsfrequenz c 

ermittelt werden. Die direkte Losung bestunde darin, die Signalfiussdiagramme in den Fig. 5 und 8 einzusetzen und 
deren Ausgangssignale nach Durchlaufen passender Glattungsfilter der Interpolationsstufe 16 (siehe Fig. 2) zuzufuh- 
ren. Das erfindungsgemasse Verfahren beinhaltet jedoch eine weitere wesentliche Vereinfachung aufgrund der Tat- 
sache, dass das FrequenzgruppenfiJter 11 lediglich eine frequenzselektive Gewichtung des Eingangssignals x vor- 

55 nimmt. Damit ist es moglich, die Eintrage in den ursprunglichen Interpolationstabellen so zu modifizieren, dass fur die 
Steuergrosse y jeweils der gleiche Wert resultiert, wenn anstelle der logarithmischen Signalleistung p des Eingangs- 
signals x die logarithmische Signalleistung q des Signais cp zusammen mit den modrfizierten Tabellen verwendet wird. 
Damit entfallt im erfindungsgemassen Verfahren die separate Berechnung der Signalleistung p, und die Verarbeitungs- 
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stufe 13 in Fig. 2 umfasst lediglich das in Fig. 8 dargestellte Signalflussdiagramm. 

Wie ebenfalls schon erwahnt wurde, weisen die aufeinanderiolgenden Signalwerte der Ausgangssignale der Ver- 
arbeitungsstufen 13 und 16 eine unerwunschte Streuung auf, die mit Glattungsfiftern 14, 15 und 17 eliminiert wird. Der 
Einsatz klassischer, linearer Tiefpassfilter ware naheliegend, ist aber wegen der damit verbundenen Verzogerungs- 

5 zeiten im erfindungsgemassen Verfahren vollig inakzeptabel. An deren Stelle wird deshalb ein nichtlineares Glattungs- 
filter gemass Fig. 9 eingesetzt, das sich nebst minimaler Verzogerungszeit noch durch wesentlich reduzierten Rechen- 
aufwand auszeichnet. Ein neuer Ausgangswert c ergibt sich dabei durch Hinzuaddieren einer Korrekturgrosse D zum 
Ausgangswert des vorhergehenden Abtastzeitpunkts. Die Korrekturgrosse D wird aus der Differenz d ermittelt, die 
sich aus dem neuen Eingangssignatwert c r und dem vorhergehenden Ausgangssignalwert ergibt. Die Grosse d wird 

10 zunachst mit einem konstanten Faktor a > 1 multipliziert. In den Glattungsfiltem 14, 15 und 17 ist der Wert von a 
beispielsweise auf 2 bzw. 3 gesetzt, und das Resultat der Multiplikation wird mit einer Sattigungsarithmetik auf den 
Wertebereich [-1, 1] begrenzt. Das Produkt w wird anschtiessend quadriert und auf einen Wert p begrenzt, und die 
Korrekturgrosse D ergibt sich durch Multiplikation des so errechneten Wertes mit der Grosse w. 

Die Wirkung des nichtlinearen Glattungsfi Iters, dessen Signalflussdiagramm in Fig. 9 gezeigt ist, wird aus der Fig. 

is 10 verstandlich, welche den Zusammenhang der intemen Grossen d und D aufzeigt. Zunachst ist zu bemerken, dass 
diese Glattungsfilter von der Normiertheit der zu filternden Signale Gebrauch machen, dass deren Wertebereich also 
das interval! [0, 1] umfasst. Damit nimmt die Differenz d Werte aus dem Intervall [-1, 1] an. Die in Fig. 10 gezeigte 
Abbildungskurve D(d) setzt sich aus funf unterschiedlichen Kurventeilen 27.1-27.5 zusammen. Fur betragsmassig 
kleine Werte der Differenz d, beispielsweise fur -0.2 < d < 0.2, hangt die Korrekturgrosse D in der dritten Potenz von 

20 der Differenz d ab; dies entspricht einem ersten Kurventeil 27.1 . Die geringfugigen Streuungen aufeinanderfotgender 
Signalwerte mit Werten aus dem Wertebereich [-0.1 ,0.1] werden damit effizient unterdruckt. Fur betragsmassig gros- 
sere Werte der Differenz d, beispielsweise fOr 0.2 < Idl < 0.5, geht die Abbildungskurve D(d) in linear verlaufende Teile 
Ober; dies entspricht einem zweiten und dritten Kurventeil 27.2 und 27.3. Bei signifikanten Veranderungen des Ein- 
gangssignals sorgen dieser Teile dafur, dass das Ausgangssignal mit nur minimaler Verzogerung folgt. Ein vierter und 

25 funfter Teil 27.4 und 27.5 der Abbildungskurve, wo eine Begrenzung auf je einen konstanten Wert erfolgt, garantiert 
schliesslich einen glatten Transitionsverlauf selbst bei extrem unstetigen Veranderungen des Eingangssignals d. 

Mit der gefilterten Schwerpunktsfrequenz c und der gefilterten Signalleistung q erfolgt in der Verarbeitungsstufe 
16 die Berechnung der Steuergrosse \\f. Wie bereits erwahnt wurde, erfolgt dieser Vorgang mittels einer zweidimen- 
sionalen Interpolation, die in Fig. 11 in einem detaillierten Schema dargestellt ist. Das Schema umfasst drei Tabellen. 

30 Die mit y 0 bezeichnete Tabelle beinhaltet die Stutzpunktwerte f Or fest vorgegebene Werte der Eingangsgrossen c und 
q. Die beiden anderen, mit dy/dc und 9\j//dq bezeichneten Tabellen enthalten die zu den Stutzpunkten passenden 
Gradientenwerte der Funktion y(c,q) in Richtung der c- und der q-Koordinaten. Der Wert der Steuergrosse y fur be- 
liebige Eingangssignalwerte c und q ergibt sich folglich naherungsweise durch 



wobei Cj und q k die zu c bzw. q nachstgelegenen Stutzpunktkoordinaten darstellen, die zugleich nicht grosser als c 
bzw. q selbst sind. Aufgrund der auf den Wertebereich [0, 1] normierten Eingangsgrossen c und q lassen sich im 
erfindungsgemassen Verfahren die Werte Cj undq k sowie (c - Cj) und (q - q k ) durch einfaches Ausmaskieren der in Fig. 
11 gezeigten Bitfelder aus den Grossen c und q ermitteln. Zur Adressierung der Tabellenwerte werden schliesslich die 
gemass Fig. 11 zusammengefassten Werte c- t und q k verwendet. 

Ein weiterer Aspekt des erfindungsgemassen Verfahrens betrifft die Verwendung optimaler Tabellenwerte in der 
zweidimensionalen Interpolation. Die Werte der Funktion ^(c.q) an den Ecken eines durch aufeinanderfolgende Stutz- 
punktkoordinaten definierten Rechtecks seien schematisch mit\|/(Cj,q k ), \j/(c k1 ,q k ), \j/(Cj,q k+1 ) und\j/(c j+1 ,q k+l ) bezeich- 
net. Dann werden im erfindungsgemassen Verfahren die Tabellenwerte 



35 



V r = ¥ 0 ( c i«Qk) + ( c " c \)< dy V /dc K,qk + (CI ■ ^k)-( d V^q)l cj qk 



(8) 



¥ 0 ( c i-q k ) = V(Cj.q k ) + [V(c k1 ,q k ) + ^(Cj.q^) - \|/(c i+v q k+1 ) - v(c jf q k )y4 



(9) 



50 



(d\j//dc)l cjqk = {[\j/(c k1 ,q k+1 ) - V(c j( q k+1 )] + [v(c j+1 ,q k ) - v(c jt q k )]}/2 



(10) 



und 



55 



(a\|//3q)l c| qk = {[^(c k1 ,q k+1 ) - v(c j+1 t q k )] + [^(c if q k+1 ) - V|/(c j( q k )]}/2 



(11) 
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verwendet. Damit werden die unvermeidlichen Interpolationsfehler gleichmassiger verteilt a!s mit den nahliegenden 
Tabellenwerten xj/(Cj,q k ), [y(c i+1 ,q k )- y(Cj,q k )J und (y(Cj,q k+1 ) - iy(Cj,q k )]. Wie ebenfalls schon erwahnt wurde, weisen die 
aufeinanderfolgenden Signalwerte y r eine unerwunschte Streuung auf , die mit dem Glattungsfilter 1 7 (vgt. Fig. 2) 
eliminiert wird. Das Ausgangssignal des Glattungsfilters 17 ist die Steuergrosse y, die in der Interpolationsstufe 5 (vgl. 

5 Fig. 1 ) zum Bestimmen von Filterparametern des Korrekturfi Iters 7 eingesetzt wird. 

Die Interpolationsstufe 5 ist im Bloc kdiag ram m der Fig. 12 detaillterter dargesteltt. Die Steuergrosse y gelangt zu 
einer Verarbeitungsstufe 18, wo aus ihr fur die nachfolgenden Interpolationen durch Ausmaskieren der in Fig. 13 
gezeigten Bitfelder eine Tabellenadresse \|/ a und eine Proportionalgrosse y f gewonnen werden. Eine Verarbeitungs- 
stufe 1 9 stellt einen 3 Bit breiten Zahter dar, dessen Zahlwert mit j bezeichnet ist. In einer Verarbeitungsstufe 20 wird 

io ein Verstarkungswert g des Korrekturfilters 7 ermittelt, und in einer Verarbeitungsstufe 21 werden Filterkoeffizienten 
kj( n > und kj(P) ermittelt. Der Zahlwert j und die interpolierten Filterparameter g, kj( n > und kj(p>, werden als Gesamtheit mit 
m bezeichnet 

Der Zahlwert j und die interpolierten Filterparameter g, kj< n > und kj<P> gelangen zum Korrekturfi Iter 7, das im Block- 
diagramm der Fig. 14 detaillierter dargesteltt ist. Es umfasst eine Verstarkerstufe 22, ein Kreuzgliedfitter 24 zur Rea- 
ls lisierung von Nullstellen und ein Kreuzgliedfitter 26 zur Realisierung von Polstellen. Der vbllstandigkeit halber sind die 
Strukturen der Kreuzgliedfilter 24 bzw. 26 in den Signalflussdiagrammen der Figuren 15 bzw. 16 ausfuhrlich wieder- 
gegeben. 

Zu jedem Abtastzeitpunkt gelangt ein interpolierter Verstarkungswert g zur Verstarkerstufe 22 (vgl. Fig. 14) und 
wird mit dem um beispielsweise 2 ms verzogerten Eingangssignal multipliziert. Die Filterkoeffizienten kj("> bzw. kj(P) 

20 gelangen zu Verarbeitungsstufen 23 bzw. 25, zu welchen auch der Zahlerwert j gef uhrt wird. Die Verarbeitungsstufen 
23 bzw. 25 sind ledigtich Schalter, welche die interpolierten Filterkoeffizienten werte entsprechend dem Zahlerwert j 
dem richtigen Filterkoeffizienten in den Kreuzgliedfiltern 24 bzw. 26 zuweisen. Den Zahlerwerten 0 bis 7 sind die Fil- 
terkoeffizienten mit den Indizes 1 bis 8 in aufsteigender Reihenfolge zugeordnet. 

Die Interpolationsstufen 20 bzw. 21 (vgl. Fig. 12) sind in den Figuren 17 bzw. 18 detailliert dargesteltt. Wie bereits 

25 erwahnt wurde, werden die aus den individuellen Lautheitsangaben ermittelten Horkorrekturdaten im erfindungsge- 
massen Verfahren als Filterparameter in einer fur die Interpolation geeigneten Formabgespeichert. Fur die Verstarkung 
handelt es sich dabei um einen logarithmischen Verstarkungswert 

30 Y = Y 0 (V a ) + VfAY(v a ) (12) 

der mit Hilfe der Tabellen y 0 und Ay, wie in Fig. 17 gezeigt, aus den Eingangsgrossen \j/ a und \j/ f interpoliert wird. Der 
Vollstandigkert halber sei erwahnt, dass im vorliegenden Falle einer eindimensionalen Interpolation die Tabelle Ay 
entfallen und der entsprechende Wert durch Differenzbildung des ausgelesenen Wertes y 0 zum nachfolgend tabellier- 
35 ten Wert jedes Mai neu berechnet werden kann. 

Zur Ermittlung des in der Verstarkerstufe 22 benotigten Verstarkungs wertes g werden aus dem Wert y anschlies- 
send durch Ausmaskieren der in Fig. 19 gezeigten Bitfelder die Adresse y a und der Proportionalwert y f gewonnen. Mit 
ihrer Hilfe wird ein Verstarkungswert 

40 

g = exp(y a ) + YfAexp(Y a ) (13) 

aus den mit exp und Aexp bezeichneten Tabellen, welche Werte der Exponentialfunktion enthalten, in einer weiteren 
Interpolation gewonnen. Die Fig. 17 stellt somit ein zweistufiges Interpolationsschema dar, das zur effizienten Ermitt- 
45 lung des benotigten Ausgangswertes wiederum von der Normiertheit der Signalwerte und darauf abgestimmter Tabel- 
len Gebrauch macht. 

Im Falle der Filterkoeffizienten werden die horbehindertenspezifischen Werte in der Form der Log-area-ratio-Ko- 
effizienten abgespeichert. Im Gegensatz zum Verstarkungswert wird zu jedem Abtastzeitpunkt jeweils nur ein Koeffi- 
zient der beiden Kreuzgliedfilter 24 und 26 neu bestimmt. Wie bereits erwahnt wurde, steuert der durch die Verarbei- 
50 tungsstufe 1 9 dargestellte Modulo-7-Zah!er den Auswahlmechanismus. Im zweistufigen Interpolationsschema der Fig. 
18 wird darum der drei Bit breite Wert des Zahlers mit der Grosse \j/ a zur aktuellen Tabellenadresse zusammengefugt. 
Fur jedes der beiden Kreuzgliedfilter 24 und 26 ergibt sich der Log -area- ratio-Koeffizient 

55 \ {v) = X 0 (v) ( ¥a ,j) + MvAX (v) ( Va ,j) (14) 

durch Interpolation mit den Tabellen Ao< v ) und AA.< V >, wobei v fur eines der Symbole n oder p steht, welche die Kreuz- 
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gliedfilter 24 und 26 zur Realisierung von Nullstellen bzw. Polstellen unterscheiden. 

Die in den Kreuzgliedfiltem 24 und 26 benotigten Filterkoeffizienten kj< n > und kj<P> werden in einer erneuten Inter- 
polation ermittelt, wobei aus jedem der Log-area-ratio- Koeffizienten X zunachst wiederum durch Ausmaskieren der in 
Fig. 20 gezeigten Bitf elder ein Adresswert X^ und eine Proportionalgrosse Xf gewonnen werden. Fur die Koeffizienten 
5 der beiden Kreuzgliedfirter 24 und 26 kann dieser Vorgang wie auch die anschliessende Interpolation selbst nachein- 
ander geschehen, was in Fig. 18 mit dem Multiplexer M angedeutet ist und insbesondere zur Folge hat, dass die mit 
tanh und Atanh bezeichneten Tabellen der Tangens-hyperbolicus-Funktion nur einmal abgespeichert sein mOssen. 
Die Filterkoeffizienten 
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k (v) = tahn(\ a (v) ) + X f (v) -Atanh(X a (v) ) (15) 



ergeben sich mit einer weiteren Interpolation, wobei fur die effiziente Implementation wiederum von der Normiertheit 
der Signalgrdssen und darauf abgestimmter Tabellen Gebrauch gemacht wird. 

15 Zusammenfassend lasst sich sagen, dass im erfindungsgemassen Verfahren zur lautheitsgesteuerten Verarbei- 

tungakustischer Signale in Schallverarbeitungsgeraten ein zu verarbeitendes akustisches Signal x vollstandig im Zeit- 
bereich verarbeitet wird. Dabei wird ausgehend vom zu verarbeitenden Signal x laufend eine fur die subjektive Laut- 
heitsempfindung Normalhorender charakteristische Steuergrosse \|/ berechnet. Das Eingangssignal x wird mit einem 
zeitabhangigen Filter 7 verarbeitet, dessen Parameter mit Hiffe der Steuergrosse \jr laufend durch Interpolation im 

20 voraus berechneter und in Tabellen abgespeicherter anwenderpezifischer Korrekturdaten neu ermittelt und auf den 
zeitabhangigen Filter 7 angewendet werden. Eine erfindungsgemasse Vorrichtung zur Durchfuhrung des Verfahrens 
weist eine Verarbeitungsstufe 4 zur iterativen Berechnung der Steuergrosse y und eine damit zeitabhangig gesteuerte 
Korrekturfilterstufe 7 aut. 

25 

Patentanspruche 

1 . Verfahren zur lautheitsgesteuerten Verarbeitung akustischer Signale in Schallverarbeitungsgeraten, dadurch ge- 
kennzeichnet, dass ein zu verarbeitendes akustisches Signal (x) vollstandig im Zeitbereich verarbeitet wird, indem 

30 ausgehend vom zu verarbeitenden Signal (x) laufend eine f 0 r die subjektive Lauthertsempf indung Normalhorender 

charakteristische Steuergrosse berechnet wird und indem das Eingangssignal (x) mit einem zeitabhangigen 
Filter (7) verarbeitet wird, dessen Parameter mit Hilfe der Steuergrosse (y) laufend durch Interpolation im voraus 
berechneter und in Tabellen abgespeicherter anwenderpezifischer Korrekturdaten neu ermittelt und auf den zeit- 
abhangigen Filter (7) angewendet werden. 

35 

2. Verfahren nach Anspruch 1, dadurch gekennzelchnet, dass das akustische Signal (x) auf eine iterative Weise 
ohne Aufteilung in Teilbandsignale verarbeitet wird. 

3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass die Steuergrosse (y) als Wurzel der auf ein 
40 begrenztes Lautheitsintervall normierten Lautheit definiert ist. 

4. Verfahren nach einem der Anspruche 1-3, dadurch gekennzeichnet, dass die Steuergrosse (\jr) laufend durch 
eine zweidimensionale Interpolation ermittelt wird, und zwar mit Hilfe zweier iterativ berechneter Grossen, von 
denen eine erste iterativ berechnete Grosse (p) ein auf einer logarithmischen Skala ausgedruckter Schatzwert fur 

45 die momentane Signalleistung und eine zweite iterativ berechnete Grosse (c) ein auf einer Barkskala ausgedruck- 

ter Schatzwert fur die Schwerpunktsfrequenz der momentanen Signalleistungsverteilung sind. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass die erste iterativ berechnete Grosse (p) mit Hilfe 
einer in eine digitale Regelschleife eingebetteten, iterativen Schatzwertberechnungseinheit erster Ordnung fur 

50 einen zeitlich exponentiell gewichteten Erwartungswert des quadrierten Eingangssignals ermittelt wird. 

6. Verfahren nach Anspruch 4 oder 5, dadurch gekennzeichnet, dass die zweite iterativ berechnete Grosse (c) 
durch Division eines iterativ ermittelten Dividenden durch einen iterativ ermitteften Divisor berechnet wird, wobei 
der Divisor ein Schatzwert fur die momentane Leistung des mit einem Frequenzgruppenfifter gewichteten Signals 

55 (cp) und der Dividend ein Schatzwert f Or die momentane Leistung des Signals (d) ist, welches daruber hinaus noch 

mit einem Barkfilter gewichtet ist, wobei die Ubertragungsfunktion des Frequenzgruppenfilters der Wurzel einer 
normierten Frequenzgruppenbreitefunktion und jene des Barkfilters der Wurzel einer normierten Tonheitsfunktion 
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entsprechen. 

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, dass sowohl der Divisor als auch der Dividend mit Hilfe 
einer in eine digitale Regetschleife eingebetteten, iterativen Schatzwertberechnungseinheit erster Ordnung fur 

5 einen zeitlich exponentieil gewichteten Erwartungswert des quadrierten Eingangssignals ermittett werden, wobei 

die Einheit zur Ermittlung des Dividenden die Regelsignale von jener des Divisors bezieht und auf ihre Signale 
anwendet. 

8. Verfahren nach Anspruch 6, dadurch gekennzeichnet, dass die Division anhand der geregelten Schatzwertgros- 
w sen berechnet und durch eine Multiplikation mit (1 - 8) angenahert wird, wobei 1 den Soltwert darstellt und I5I < 

1 gilt. 

9. Verfahren nach einem der Anspruche 5-8, dadurch gekennzeichnet, dass die zur Regelung der iterativen Schatz- 
wertberechnungseinheit benotigten Skalierungsgrdssen sowie die zur Nachfuhrung des logarithm ischen Schatz- 

is wertes benotigten inkrementellen Ande rungs werte aus im voraus abgespeicherten Tabellen (S, A, Ap) ausgelesen 

werden. 

10. Verfahren nach Anspruch 9, dadurch gekennzeichnet, dass das Auslesen aus auf solche Weise organisierten 
Tabellen erfolgt, dass sich die Tabellenindizes zum Auffinden der gesuchten Grossen durch einfaches Ausmas- 

20 kieren von Bitfeldern aus der noch ungeregetten Schatzwertgrosse (v) und der logarithm ischen Schatzwertgrosse 

(p) gewinnen lassen. 

11. Verfahren nach einem der Anspruche 1-3, dadurch gekennzeichnet, dass die Steuergrosse laufend durch 
eine zweidimensionale Interpolation ermittett wird, und zwar mit Hilfe zweier iterativ berechneter Grossen, von 

25 denen eine erste iterativ berechnete Grosse (q) ein auf einer logarithmischen Skala ausgedrOckter Schatzwert fur 

die momentane Leistung eines mit einem Frequenzgruppenfilter gewichteten Signals (<p) ist, wobei die Gewichtung 
durch Anderung der Eintrage in der ursprunglichen Interpolationstabelle kompensiert wird, und eine zweite iterativ 
berechnete Grosse (c) ein auf einer Barkskala ausgedrOckter Schatzwert fur die Schwerpunktsfrequenz der mo- 
mentanen Signalleistungsverteilung ist. 

30 

12. Verfahren nach einem der Anspruche 4-11, dadurch gekennzeichnet, dass die Steuergrosse und/oder die 
erste iterativ berechnete Grosse (p bzw. q) und/oder die zweite iterativ berechnete Grosse (c) mit einem nichtli- 
nearen Filter geglattet werden, und zwar derart, dass sich ein neuer Ausgangswert durch Addition eines Korrek- 
turwertes (D) zum vorherigen Ausgangswert ergibt, dass sich dieser Korrekturwert (D) aus der Differenz (d) zwi- 

35 schen neuem Eingangssignal und vorherigem Ausgangssignal berechnet und dass der Korrekturwert (D) fur kleine 

Betragswerte (Idl) der Differenz (d) in der dritten Potenz von dieser Differenz (d) abhangt, fur mittlere Betragswerte 
(Idl) der Differenz (d) linear von dieser Differenz (d) abhangt und fur grosse Betragswerte (Idl) der Differenz (d) 
konstant ist. 

40 13. Verfahren nach einem der Anspruche 4-12, dadurch gekennzeichnet, dass die Interpolation der Steuergrosse 
(y) mit derart organisierten Tabellen erfolgt, dass sowohl der Tabellenindex zum Auffinden des Stutzpunktwerts 
und der inkrementellen Zuwachsgrossen in beiden Dimensionen als auch die Proportionalgrossen, mit denen die 
inkrementellen Zuwachsgrossen vor der Addition zum Stutzpunktwert multipliziert werden, durch einfaches Aus- 
maskieren von Bitfeldern aus den iterativ berechneten Grossen (p bzw. q; c) gewinnbar sind. 

45 

1 4. Verfahren nach Anspruch 1 3, dadurch gekennzeichnet, dass in den Tabellen zur zweidimensionalen I nterpolation 
der Steuergrosse (y) optimierte Werte gemass den Formeln 

so V 0 ( c i^k) = V(Ci>q k ) + [V( c j + i.q k ) + V( c j.q k+ i) * V(c kl .q k+1 ) - V(Cj,q k )]/4 (9) 

(dy/dc)^ qk = «v(c. +1 ,q k+1 ) - y(Ci,q k+1 )] + [ V (c kl ,q k ) - y{c v q k )]}/2 (10) 
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(ay/dq)l dqk = {[V(c i+1 .q k+1 ) - v(c w .q k )] + [v(c,,q k+1 ) - y{c v q k )]}/2 (11) 

verwendet werden. 

5 

15. Verfahren nach einem der Anspruche 1-14, dadurch gekennzeichnet, dass die zur Interpolation der anwender- 
spezifischen Korrekturdaten in Tabellen gespeicherten Werte ats Verstarkungswerte im logarithmischen Bereich 
und als Filterkoeffizienten im Log-Area- Ratio- Bereich abgelegt sind. 

10 16. Verfahren nach Anspruch 15, dadurch gekennzeichnet, dass die Interpolation der anwenderspezifischen Kor- 
rekturdaten mit derart organisierten Tabellen erfolgt, dass der Tabellenindex zum Auffinden des Stutzpunktwertes 
sowie der Tabellenindex zum Auffinden der Proportionalgrosse, mit der die Differenz zwischen dem nachstfolgen- 
den Stutzpunktwert und dem Stutzpunktwert selbst vor der Addition zum Stutzpunktwert multipliziert wird, durch 
einfaches Ausmaskieren von Bitfeldern aus der Steuergrosse (y) gewonnen werden. 

15 

17. Verfahren nach Anspruch 15 Oder 16, dadurch gekennzeichnet, dass der Verstarkungswert aus dem interpolier- 
ten logarithmischen Verstarkungswert und die Filterkoeffizienten aus den interpolierten Log -Area- Ratio-Koeffizi- 
enten wieder durch Interpolation mit abgespeicherten Tabellen der Exponentialfunktion und der Tangenshyperbo- 
likusfunktion sowie Tabellen der inkrementellen Zuwachsgrossen dieser Funktionen ermittelt werden. 

20 

18. Verfahren nach Anspruch 17, dadurch gekennzeichnet, dass die Interpolation mit derart organisierten Tabellen 
erfolgt, dass sich die Tabellenindizes zum Auffinden der Stutzpunktwerte und der inkrementellen Zuwachsgrossen 
sowie die Proportionalgrossen, mit denen die inkrementellen Zuwachsgrossen vor der Addition zu den StOtzpunkt- 
werten multipliziert werden, durch einfaches Ausmaskieren von Bitfeldern des interpolierten Verstarkungswertes 

25 und der interpolierten Log-Area-Ratio-Koeffizienten gewonnen werden. 

19. Verfahren nach einem der AnsprOche 15-18, dadurch gekennzeichnet, dass der Verstarkungswert in jedem Abt- 
astintervall und von den Filterkoeffizienten in jedem Abtast interval I jeweils nur die Koeffizienten eines Pol-/Nullst- 
ellenpaares neu ermittelt werden, wobei fur die Erneuerung der Filterkoeffizienten eine feste, gleichmassige Rei- 

30 henfolge gilt. 

20. Verfahren nach einem der Anspruche 1 -1 9, dadurch gekennzeichnet, dass das Eingangssignal zum vorgenann- 
ten zeitabhangigen Filter derart verzogert wird, dass die uber die Berechnung der vorgenannten Grosse stets 
neu ermittelten Filterkoeffizienten und Verstarkungswerte zeitgerecht auf das der Berechnung zugrunde liegende 

35 Signal angewendet werden. 

21. Vorrichtung zur Durchfuhrung des Verfahrens nach Anspruch 1, gekennzeichnet durch eine Verarbeitungsstufe 
(4) zur iterativen Berechnung der Steuergrosse (y) und eine damit zeitabhangig gesteuerte Korrekturfilterstufe (7). 

40 22. Vorrichtung nach Anspruch 21 , gekennzeichnet durch eine zweidimensionale Interpolationsstufe (16) zur Ermitt- 
lung der Steuergrosse (\|/) aus einer Signalleistung (q) und einer Schwerpunktsfrequenz (c). 

23. Vorrichtung nach Ansprch 21 oder 22, gekennzeichnet durch ein Frequenzgruppenfilter (11) und ein Barkfilter 
(12) zur Ermittlung von gefilterten Signalen (<p, #) aus einem Eingangssignal (x). 

45 

24. Vorrichtung nach Anspruch 23, dadurch gekennzeichnet, dass das Frequenzgruppenfilter und das vorgenannte 
Barkfilter als Rekursiv-filter ausgelegt sind. 

25. Vorrichtung nach einem der Anspruche 21-24, gekennzeichnet durch eine Schatzwertberechnungseinheit (13) 
50 zur Berechnung der Signalleistung (q) und der Schwerpunktsfrequenz (c) aus den gefilterten Eingangssignalen 

(<P. 

26. Vorrichtung nach einem der Anspruche 21 -25, gekennzeichnet durch Glattungsfifter (14, 15, 17) zur Eliminierung 
unerwunschter Streuung aufeinanderfolgender Signaiwerte (c r , q r , y r ). 



55 



27. Vorrichtung nach einem der Anspruche 21 -26, gekennzeichnet durch eine Serieschaltung einer Verstarkungs- 
stufe (22), einer Kreuzgliedfitterstufe (24) zur Realisierung von Nullstellen und einer Kreuzgliedfilterstufe (26) zur 
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Realisierung von Polstellen. 

28. Vorrichtung nach einem der AnsprOche 21-27, gekennzeichnet durch zweistufige Interpolationsstufen zur Er- 
mittlung des Verstarkungswertes (g) und der Koeffizienten (k|<°) und kj<P)) des Korrekurfilters (7) aus der Steuer- 
grosse 

29. Vorrichtung nach einem der AnsprOche 21-28, gekennzeichnet durch eine Signalverzogerungseinheit (6) zur 
Synchronisation des Eingangssignals (x) hinsichtlich der Verarbeitung mil dem Korrekturfifter (7), dessen Filter- 
parameter aus dem Eingangssignal (x) abgeleitet werden. 
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